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ABSTRACT 

Through a series of observations with the Austraha Telescope Compact Array 
we have monitored the variabihty of ground-state hydroxyl maser emission from 
G12. 889+0. 489 in all four Stokes polarisation products. These observations were mo- 
tivated by the known periodicity in the associated 6.7-GHz methanol maser emission. 
A total of 27 epochs of observations were made over 16 months. No emission was seen 
from either the 1612 or 1720 MHz satellite line transitions (to a typical five sigma up- 
per limit of 0.2 Jy). The peak flux densities of the 1665 and 1667 MHz emission were 
observed to vary at a level of ^^20% (with the exception of one epoch which dropped 
by ^40%). There was no distinct flaring activity at any epoch, but there was a weak 
indication of periodic variability, with a period and phase of minimum emission sim- 
ilar to that of methanol. There is no significant variation in the polarised properties 
of the hydroxyl, with Stokes Q and U fiux densities varying in accord with the Stokes 
I intensity (linear polarisation, P, varying by ^20%) and the right and left circularly 
polarised components varying by ^33% at 1665-MHz and ^38% at 1667-MHz. These 
observations are the first monitoring observations of the hydroxyl maser emission from 
G12.889-H0.489. 
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1 INTRODUCTION 

A site of hydroxyl maser emission, G12. 889-1-0. 489, initially 
found in a search towards IRAS 18089-1732 (|Cohen et al.1 
Il988l ). wa s found to also be a prominent methanol maser at 
6.7 GHz (lMentenlll991^. V a riabili ty of the methanol was 



identified by ICaswell et aL ll 19951) and l ed to the exten- 



sive monitoring o f iGoedhart et all (l2004f) and the discov- 



ery of periodicity (|Goedhart et al.l 20091 ). 6.7-GHz methanol 
masers exclusively trace the formation of high-mass stars 
and are hence closely studied to gain insight into the largely 
unknown processes of high mass star formation. The discov- 
ery of periodicity in the methanol maser emission indicates a 
periodic variation in the inputs of the maser emission process 
which has significant implications both for the mechanism 
of maser emission and the nature of high-mass star forma- 
tion, indicating periodic processes such as the interaction of 
winds from binary pre-m ain sequence high mass stars (e.g. 
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Goedhart et al.l ( 20091 ) found intensity variations in the 



methanol emission with a period of nearly 30 days that ap- 
peared stable, although variations from cycle to cycle in the 
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peak amplitude of the fiare were apparent. The flaring fea- 
tures peaked anywhere within an 11-day window, but the 
phase of the minima was stable. An amplitude variation with 
the same period is present in features at different velocities 
and at both 6.7 GHz and 12.2-GHz transitions. Delays of up 
to six days were found between individual 6.7-GHz features 
and a one day delay was found between the 12.2-GHz and 
6.7-GHz flares at the same velocities. 

Hydro xyl observations since the discovery of the maser 
bv lCoheiTe t al. (1988) have i ncluded the positio ning obser- 
vations of PCasweU (|1998D a nd I Argon et al.1 (|2000|) . and a de- 
tailed polarization study bv lSzvmczak fc GerardI (|2009l ). Ap- 
parent variations in total intensity exceed 20%, but intervals 
between observations were typically longer than one year 
and precise comparison is not possible with the significantly 
different instrumentation. Observations with the Parkes ra- 
dio telescope (Caswell et al. in prep.) measured the hydroxyl 
emission in 2004 and 2005, finding peak flux densities of 8 
Jy at 1665 MHz and 1.8 Jy at 1667 MHz. The Parkes spec- 
tra are in agreement at the tw o epochs to within 10% and 
similar values were obtained bv lSzvmczak fc GerardI (|2009| ) 
in observations made in 2003. Although the amplitudes have 
shown some variation, the velocities of the features are con- 
sistent over all previous observations. 
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The majority (~80%) of hydroxyl masers in regions of 
high-mas s star formati on have associated methanol maser 
emission ()Caswelllll998l l. with many sites exhibiting a close 
spatial coincidence, implying a common masing gas and 
pumping source. Models of maser pumping suggest that 
both species are pumped by infrared emission from dust 
surrounding the h i gh mass pre-main s equence obj ect (e.g. 
iMoore et all Il988l : ICragg et al.l l2005l : iGravl I2OO7I ). Hence 
there is an expectation that any variability in emission of 
one species will correlate with variability in the other. Early 
evidence for hydroxyl maser variability fe.g. iRobinson et al.l 
I1970D w as followed by several more focused variability stud - 
ies (e.g. ISullivan fc Kersthol'3ll976l : IClegg fc Cordeslll99lh . 
but these showed no clear periodicity of hydroxyl masers 
in star forming regions. More recent investigations of possi- 
ble common flaring behaviour between the methanol and 
the ground-state hydroxyl ma sers has been inconclusive 
fe.g. iMacLeod fc Gavlardlll996l ). but co rrelated variability 
of ex cited-state hydroxyl with methanol (lAl-Marzouk et al] 
l2012f ) and formaldehyde with methanol (|Arava et al. 2010l ) 
has been found recently. With methanol emission from 
G12. 889-1-0. 489 established to be periodic in nature, we were 
motivated to initiate monitoring observations of the associ- 
ated hydroxyl emission. 



2 OBSERVATIONS 

Observations were made with the Australia Telescope Com- 
pact Array (ATCA), commencing with a full synthesis image 
on 2010 May 20 consisting of four cuts of 10 minutes spread 
over eight hours. This was followed by 26 epochs over 1.5 
years with typical integration times of 40 minutes on source. 
These observations were divided into two groups, the first 
with one to two observations per week for ~10 weeks, the 
second with one observation approximately every 24 hours 
for seven days. Observations were centred at 1612.2310 MHz, 
1665.4018, 1667.3590 MHz and 1720.5300 MHz, corrected 
for the motion of the local standard of rest. The ATCA cali- 
brator 1830—210 was used for phase calibration for the first 
group of observations and ATCA calibrator 1829—207 for 
the second group (with positional uncertainties of <0.15" 
and <0.0l" respectively). Both groups used the primary 
ATCA flux density cahbrator PKS 1934-638, with flux den- 
sities bootstrapped to this (with a resultant systematic un- 
certainty of <2%). All observations wer e made with the new 
Compact Array Broad Band backend (| Wilson et al.|[201ll ). 
adopting the CFB lM-0.5k mode (2048 channels over 1 MHz 
giving 0.5 KHz channel spacings), obtaining all four polarisa- 
tion products. The velocity channel separations were 0.091, 
0.088, 0.088 and 0.085 km s"^ for the 1612, 1665, 1667 and 
1720 MHz transitions respectively. The full listing of the ob- 
servations and noise levels is given in Table[T] Antenna 4 was 
not used for epochs seven to 18 due to the upgrade of the 
low frequency (1—3 GHz) receiver. Observations were made 
within the LST range 14:30 to 22:00 hrs and predominantly 
in the 6-km array conflgurations (6B and 6C), but eight 
epochs were observed with compact configurations (see Ta- 
ble[T]for details). There were no systematic variations in flux 
density for the observations with the short baseline (com- 
pact) array configurations compared with the long baseline 
(6-km) configurations. The data were reduced and processed 
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Figure 1. Range in flux density variation in Stokes I for 
G12. 8894-0. 489 at each ground-state hydroxyl transition. The 
solid line is the averaged spectrum, constructed from the me- 
dian value from all epochs for each spectral channel. The two 
dashed lines show the net extreme values (having excluded the 
most extreme noise- biased values). 



with the MIRIAD softwa re package using standard techniques 
(|Sault fc Killeenll2004l '). Following the positioning observa- 
tion, spectral profiles were obtained with the task uvspec 
and a time series of peak flux density of maser features was 
determined through fitting a Gaussian in the spectral do- 
main with minimisation for each epoch. 



3 RESULTS 

The position of strongest emission at 1665 and 1667 
MHz from the initial (full synthesis imaging) observation 
was found to be RA(J2000) 18''11'"51.45" and Dec(J2000) 
-17°31'29.7" with a positional error of sSO.4". This is coin- 
cident (to within the err ors) with an earlier ATCA mea- 
surement (jCaswelll Il998h an d with a Very La rge Array 
measurement at 1665 MHz (|Argon et all I2OOOD . Position 
estimates of the strongest methanol fe atures at 6.7 and 
12.2 GHz (|Caswellll2009l : IXu et al.ll201ll ) are 18'ni'"51.40=, 
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Table 1. Details of observations. ' - ' mark observations where data were unusable. The first observation, denoted 'imaging' was a 
synthesis imaging observation as discussed in the text. 



epoch 


Gregorian 
Date 


Modified 
Julian 
Date 


Array 
Config. 


Integration 
Time 
(mins) 


1612 
(mjy) 


<j 

1665 
(mJy) 


rms 

1667 
(mJy) 


1720 
(mJy) 


Li I LyJiSjLl L\j 


2(11 n-n'i-2n 


55336, 


,5 


6C 


32 


57 


64 


49 


64 


\ 


901 n 94 

J.VJ L/tJ .^rt 


55340, 


,5 


6C 


40 


44 


59 


49 


56 


2 


2(11 (i-ftfi-ns 


55363, 


,5 


6C 


40 


34 


35 


42 


41 


3 


901 nfi 07 


55367, 


,5 


6C 


40 


69 


42 


39 


58 


4 


901 O-Ofi-I 

J.VJ L/U J-VJ 


55370, 


,5 


6C 


40 


32 


41 


34 


49 


5 


901 Ofi 1 4 

^VJ J.VJ L/U -Lrt 


55374, 


,5 


6C 


40 


34 


34 


34 


44 


Q 


201 (l-flfi-l S 


55378, 


,5 


6C 


70 


29 


32 


33 


41 




901 Ofi 97 


55387, 


,5 


6C 


60 


38 


40 


49 


42 


g 


901 0-nfi-9Q 


55389, 


,5 


6C 


40 


50 


51 


57 


65 


q 


901 07 09 


55392, 


,5 


6C 


40 




55 


49 




10 


901 0-07-04 


55394, 


,5 


6C 


40 


45 


46 


53 


51 


w 


901 07 0^^ 


55398, 


,5 


1.5D 


45 


40 


62 


62 


45 


12 


901 0-07-1 ^ 


55405, 


,5 


F,W352 

1 ' V V rJlJZj 


40 


200 


137 


141 


259 


13 


2010-07-18 


55408, 


,5 


EW352 


40 


184 


219 


131 


254 


14 


2010-07-21 


55411, 


,5 


EW352 


45 


113 


90 


91 


102 


15 


2010-07-27 


55417, 


,5 


H168 


80 


52 


53 


57 


57 


16 


2010-07-29 


55419, 


,5 


H168 


45 


56 


61 




68 


17 


2010-08-01 


55422, 


,5 


H168 


40 


82 


90 


91 


96 


18 


2010-08-05 


55426, 


,5 


H168 


72 


46 


66 


55 




19 


2011-08-26 


55812, 


,5 


6B 


40 


33 


37 


36 


36 


20 


2011-08-27 


55813, 


,5 


6B 


40 


42 


37 


36 


38 


21 


2011-08-29 


55815, 


,5 


6B 


40 


34 


39 


39 


35 


22 


2011-08-30 


55816, 


,5 


6B 


40 


37 


41 


35 


35 


23 


2011-08-31 


55817, 


,5 


6B 


40 


33 


35 


37 


35 


24 


2011-09-01 


55818, 


,5 


6B 


40 


48 


39 


34 


36 


25 


2011-09-02 


55819, 


,5 


6B 


40 


43 


38 


43 


39 


26 


2011-09-03 


55820, 


,5 


6B 


40 


31 


37 


34 


35 



-17°31'29.6" and 18''11™51.396=, -17°3l'29.9l" respec- 
tively, a nominal ofltset from the hydro xyl of 0.7"±0.6 ". At 
an astrometric distance of 2.3±0.1 kpc l|Xu et al.ll201ll '). this 
small angular offset would correspond to a physical offset of 
1500 AU (~10 light days). We found that the peak flux den- 
sities of the main features in both the 1665-MHz and 1667- 
MHz transitions showed variability, but by less than 20±3% 
from the maximum peak flux density with the exception 
of one epoch in the 1665-MHz transition and one epoch in 
the 1667-MHz transition, where emission dropped by 31±3% 
and 38±3% respectively. The averaged spectra and extrema 
for the ground state transitions are given in Figure [T] The 
average spectrum is the combination of the median value at 
each spectral channel taken from the values of all epochs 
of observations. Similarly the extrema are the highest and 
lowest values for a spectral channel across all epochs, exclud- 
ing the furthest outliers (to avoid bias by any extreme noise 
fluctuations). The Gaussian fits used to construct the time 
series had errors of 0.01 kms~^ for both of the 1665-MHz fea- 
tures (found to be centred at 33.01 kms^^ and 35.26kms~^) 
and <0.01kms"\ 0.02kms~\ 0.02 kms"^ and 0.03 kms~^ 
for the four (weaker) 1667 MHz features (at 33.46 km s~^, 
32.82kms~^, 34.79kms~^ and 31.53kms~^ respectively). 
Any variation in the fitted peak velocity was within the er- 
rors and no systematic shift in the peak velocity of features 
was found. There was also no consistent shift in velocity for 
the epochs which had minima in the flux density. The time 
series for the main features at each frequency are shown in 
Figures [2] and m No emission was detected from the satellite 



line transitions at any of the epochs (to ~0.2 Jy 5(t limit). 
Additionally, for the last eight epochs the excited state tran- 
sitions at 6030 MHz and 6035 MHz were observed, but also 
did not have any detected emission (to ~0.2 Jy ba limit). 

A Lomb-Scarele per iodogram (following the algorithm 
of|P ress fc Rvbickil [1989!) was used to search for any period- 
icity in each of the six features shown in the time series, with 
results displayed in Figure[4l At low significance there is a 
suggestion of a 0.03 to 0.04 cycles per day peak which would 
correspond to a period of 25 to 30 days, coincident with the 
period of the 6.7-GHz methanol maser. It is most noticeable 
in the 33.0 and 35.3 kms~^ peak features of the 1665-MHz 
transition which have power peaks at 0.034 cycles per day 
(29.4 days). This power peak was not affected by excluding 
the two highest flux density epochs from the periodogram. 
In contrast, after exclusion of the two lowest flux density 
epochs, the power peak was less evident for the 33kms~^ 
feature and absent for the other features. This implies, com- 
parable to the methanol, that the suggested periodicity is 
predominantly in the minima of the emission. 

The features in the 1665-MHz Stokes I spectrum at 
31.6 km s~^, 33.0 km s~^ and 34.0 km s~^ were found on av- 
erage to be 91%, 55% and 70% linearly polarised respec- 
tively. The features in the 1667-MHz Stokes / spectrum at 
31.5kms"\ 32.9kms~\ 33.4 km s"^ and 34.8 km s"^ were 
found on average to be 84%, 72%, 23% and 2% linearly 
polarised respectively. The averaged and extreme spectra of 
the linear polarisation are shown in Figure[5l The right hand 
circular polarisation (RHCP) and left hand circular polari- 
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Figure 2. Time series for the two brigiitest features at 1665 MHz and the four brightest features at 1667 MHz in Stokes / for the 
first group of observations, spanning 90 days. Errors are combination of the rms noise and the errors in the Gaussian fit in the spectral 
domain. Solid fdashed ) grey vertical lines show epochs of extrapolated minima (and half a period offsets) of methanol emission from 
iGoedhart et al.i l|2009l V 



sation (LHCP) features, derived from RHCP = (/ + V)/2 
and LHCP = (I — V )/2, are shown in Figure[6] Spectra from 
lArgon et al] l|2000l ) taken in 1993 January, limited to 1665 
MHz, and only RHCP and LHCP, show the same features as 
our spectra but typically at half the intensity. Archival un- 
published data from Parkes confirm the lower intensity near 
that epoch (1993 July). The Argon et al. VLA data sug- 
gest that some spectral features are spatial blends of several 
nearby components (although not well resolved by the VLA 
beam), but the separations of all major features in the ve- 
locity range 32 to 36kms~^ are less than 0.15". 

In our monitoring sessions, of all the circularly polarised 
features, the RHCP features were found to vary by up to 
33±3% from the maximum and the LHCP features by up to 
38±3% from the maximum. Our measurement of the feature 
found at 31.6 kms~^ in both the 1665 and 1667 MHz tran- 
sitions is of particular note. It is an outlying feature both 
in velocity and spatially. The 1665 MHz emission has an 
oflFset position of RA(J2000) 18''11™51.53" and Dec(J2000) 
-17°31'27.3" (the 1667 MHz RA(J2000) 18''11'"51.54= and 
Dec(J2000) -17°31'28.6"). Both transitions have high frac- 
tional linear polarisation. The polarisation position angle of 
this feature at 1665-MHz varies between 6° and 16° with a 
median angle of 12°. The averaged spectrum has a position 
angle of 12°. The position angle of this feature at 1667-MHz 
varies between 0° and 10° with a median angle of 5°. The 
averaged spectrum has a position angle of 5°. The error in 
the position angle due to noise is estimated to be 5°. 



4 DISCUSSION 

As stated in the introduction, the methanol maser counter- 
part of G12. 889+0. 4 89 has quasi-periodic f laring at both 6.7 
GHz and 12.2 GHz (|Goedhart et all 120091 '). There are time 
delays between the two frequencies and the flaring period 
occurs anywhere within an 11 day window, but minima are 
regular (with a period of 29.5 days). This periodic variabil- 
ity was ascribed by the authors as likely being due either 
to variation in the radio continuum emission (the seed pho- 
tons for the maser emission) or variation in the infrared 
emission (pumping the maser). We note that, extrapolating 
from the work of Gocdhart e t al. (2009,) , our range of epochs 
should have covered flaring activity and minima (three min- 
ima should occur in the first session of observations and 
the second session should be wholly between two adjacent 
minima). Goedhart et al. show that the methanol exhibits 
flux density variations of 60 to 70%. With the exception of 
one epoch, the hydroxyl emission we measure varies by only 
20%. The current data have no strong indication of peri- 
odicity, but do suggest at low significance variation with a 
comparable period to the methanol, with the minimum at 
MJD 55379 nominally agreeing with the extrapolated min- 
ima epoch of the methanol. Further observing epochs will 
be required to make firm conclusions on the periodicity of 
variation. 

The polarisati on properties we obs e rve a re compara- 
ble to the work of ISzvmczak fc GerardI (|2009l ), who mea- 
sured the full polarisation properties at 1665 and 1667 MHz 
with all four Stokes parameters with the Nancay telescope in 
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Figure 3. Time series for tiie two brightest features at 1665 MHz 
and the four briglitest features at 1667 MHz in Stokes / for the 
second group of observations, spanning 8 days. Errors are combi- 
nation of the rms noise and the errors in the Gaussian fit in the 
spectral domain. Solid grey vertical line shows an epoch of the 
extrapolated minima of methanol emission from ICoedhart et al.l 
J200g ). 



2003. They found linear polarization of 87.3% and 76.7% in 
two 1665-MHz features at 31.59 (the kinematic and spatial 
outlier mentioned previously) and 32.76 kms~^ respectively, 
and both features had low circular polarisation at the epoch 
of the observations (14.6% and 8.9% respectively). If the per- 
sistently strong features, 1665-MHz RHCP at 35.3 km s^^ 
and the 1665-MHz LHCP at 33.1kms^^, are a Zeeman pair, 
the implied magnetic field strength is +3.7 mG. This is cor- 
roborated with the 1667-MHz emission where the brightest 
RHCP at 34.8 km s"^ and the brightest LHCP feature at 
33.4 km s~^, also give an implied magnetic field strength of 
+3.7 mC. 

A qualitative explanation for some of the methanol 
meiser sources with long periods (e.g. 9.62-1-0.20 and 
188.95-1-0.89) is that of a colliding wind bin ary system 
l|van der Walt et all l2009l : Ivan der Wal^ 120111 ). a system 
which could periodically alter either the seed flux or pump- 
ing mechanism of the maser emission (although for the ex- 
ample sources in the work of van der Walt et al. the seed flux 
is the more likely in view of their long periods). The hydro- 
gen winds of the binary pre-main sequence stars either heat 
the circumstellar dust or cause additional ionisation of hy- 
drogen surrounding the forming high mass star. In the case 
of G12. 889-1-0. 489, an upper limit to continuum emission of 
1 mjy has been established (Walsh e t al. 1998), but if there 
is a very weak region of ionised hydrogen providing the seed 
photons, the possible offset location of the hydroxyl relative 
to the methanol (by ~1500 AU) may place the hydroxyl suf- 
ficiently far away so as to be minimally affected, explaining 



the lower prominence of periodicity in the hydroxyl emission. 
On the other hand the optical depth at the lower frequency 
of the hydroxyl emission could be such that the seed radia- 
tion does not respond as it does at the methanol frequency. 
An aspect of this model is that we may expect to see ro- 
tation of polarisation angle of features with the passing of 
shocks associated with the colliding winds (and the flares 
of maser emission). As mentioned in the results we have 
one feature at 31.6 km s~^ with high linear polarisation, but 
this does not show significant variation in polarisation angle 
(variation is within the noise). 

An alternative theory to a ccount for periodic ity in 
masers has been put forward by lArava et al.l (|2010l ) . This 
was proposed to explain the periodicity in formaldehyde and 
methanol maser emission in IRAS 185664-0408 through cir- 
cumbinary disk accretion. In this model the accreting ma- 
terial heats the dust, increasing the photons pumping the 
maser. In this variety of model, an offset location could also 
be sufficient to diminish the effect on the hydroxyl emission. 

An important aspect of G12. 889-1-0. 489 is the shortness 
of the periodicity (29.5 days), which any model for the sys- 
tem must account for. Additionally, the suggestion that the 
minima of hydroxyl emission may coincide with the minima 
of the methanol implies that the periodicity may arise from 
a quenching or suppression mechanism, rather than flaring. 
Continued monitoring of this maser source across the vari- 
ous maser transitions is required to provide further insight. 



5 SUMMARY 

We have obtained the first monitoring observations of 
G12.889-I-0.489 with the Austraha Telescope Compact Ar- 
ray, finding variability at the level of 20% on timescales 
of days. Despite no distinct fiaring characteristics, the hy- 
droxyl variability showed possible periodicity similar to that 
of methanol, with the minimum of emission at the same 
phase. It will require a full programme of sensitive daily 
monitoring of G12. 889-1-0. 489 at the ground-state hydroxyl 
transitions to fully test the apparent periodic behaviour in 
the variability, and thus assess whether the mechanism caus- 
ing periodicity in the methanol maser counterpart is shared 
with the hydroxyl. 
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Figure 4. Results of Lomb-Scargle periodograms for the two 1665 
MHz features and four 1667 MHz features. The power spectra 
were oversampled by a factor of five and have been smoothed with 
a bezier function for clarity. There is a suggestion of a peak at 
0.0,3 to 0.04 cycles per day, which would correspond to a period of 
25 to 30 days, coincident with the period of the 6.7-GHz methanol 
maser. 
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Figure 5. Range in flux density variation in Stokes Q and U 
and Linear Polarisation (P) for G12.889+0.489 at the 1665 and 
1667 MHz hydroxyl transitions. The solid line is the averaged 
spectrum, constructed from the median value from all epochs for 
each spectral channel. The two dashed lines show the net extreme 
values (having excluded the most extreme noise-biased values). 
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Figure 6. Range in flux density variation in RHCP and LHCP for 
G12.889-I-0.489 at the 1665 and 1667 MHz hydroxyl transitions. 
The solid line is the averaged spectrum, constructed from the 
median value from all epochs for each spectral channel. The two 
dashed lines show the net extreme values (having excluded the 
most extreme noise- biased values). 



